非全局Lipschitz条件下跳适应向后Euler方法的强收敛性分析

本文研究跳适应向后Euler方法求解跳扩散随机微分方程在非全局Lipschitz条件下的强收敛性.通过克服方程非全局Lipschitz系数给收敛性分析带来的主要困难,我们成功地建立了跳适应后向Euler方法的强收敛性结果并得到相应的收敛率.最后,我们通过数值试验对前文所得理论结果做进一步的验证.     

几乎各向同性的高维空间分数阶扩散方程的分块快速正则Hermite分裂预处理方法

一类空间分数阶扩散方程经过有限差分离散后所得到的离散线性方程组的系数矩阵是两个对角矩阵与Toeplitz型矩阵的乘积之和.在本文中,对于几乎各向同性的二维或三维空间分数阶扩散方程的离散线性方程组,采用预处理Krylov子空间迭代方法,我们利用其系数矩阵的特殊结构和具体性质构造了一类分块快速正则Hermite分裂预处理子.通过理论分析,我们证明了所对应的预处理矩阵的特征值大部分都聚集于1的附近.数值实验也表明,这类分块快速正则Hermite分裂预处理子可以明显地加快广义极小残量(GMRES)方法和稳定化的双共轭梯度(BiCGSTAB)方法等Krylov子空间迭代方法的收敛速度.     

曲边区域上的多边形网格间断有限元离散及其多重网格算法

本文利用多边形网格上的间断有限元方法离散二阶椭圆方程,在曲边区域上,采用多条直短边逼近曲边的以直代曲的策略,实现了高阶元在能量范数下的最优收敛.本文还将这一方法用于带曲边界面问题的求解,同样得到高阶元的最优收敛.此外我们还设计并分析了这一方法的\linebreakW-cycle和Variable V-cycle多重网格预条件方法,证明当光滑次数足够多时,多重网格预条件算法一致收敛.最后给出了数值算例,证实该算法的可行性并验证了理论分析的结果.     

Elliptic Reconstruction and a Posteriori Error Estimates for Fully Discrete Semilinear Parabolic Optimal Control Problems

This article studies a posteriori error analysis of fully discrete finite element approximations for semilinear parabolic optimal control problems. Based on elliptic reconstruction approach introduced earlier by Makridakis and Nochetto [25], a residual based a posteriori error estimators for the state, co-state and control variables are derived. The space discretization of the state and co-state variables is done by using the piecewise linear and continuous finite elements, whereas the piecewise constant functions are employed for the control variable. The temporal discretization is based on the backward Euler method. We derive a posteriori error estimates for the state, co-state and control variables in the $L^\infty(0,T;L^2(\Omega))$-norm. Finally, a numerical experiment is performed to illustrate the performance of the derived estimators.     

An $L^{\infty}$ Second Order Cartesian Method for 3D Anisotropic Interface Problems

A second order accurate method in the infinity norm is proposed for general three dimensional anisotropic elliptic interface problems in which the solution and its derivatives, the coefficients, and source terms all can have finite jumps across one or several arbitrary smooth interfaces. The method is based on the 2D finite element-finite difference (FE-FD) method but with substantial differences in method derivation, implementation, and convergence analysis. One of challenges is to derive 3D interface relations since there is no invariance anymore under coordinate system transforms for the partial differential equations and the jump conditions. A finite element discretization whose coefficient matrix is a symmetric semi-positive definite is used away from the interface; and the maximum preserving finite difference discretization whose coefficient matrix part is an M-matrix is constructed at irregular elements where the interface cuts through. We aim to get a sharp interface method that can have second order accuracy in the point-wise norm. We show the convergence analysis by splitting errors into several parts. Nontrivial numerical examples are presented to confirm the convergence analysis.     

A Space-Time Polynomial Collocation Method for Solving 3D Burgers Equations北大核心CSCD

Burgers方程是一类应用广泛的非线性偏微分方程,方程中的非线性项难以处理。该文提出一种新的时空多项式配点法——多项式特解法求解三维Burgers方程。求解过程分为两步:第一步,对三维Burgers方程中的线性导数项(包括时间导数项),求出相应的多项式特解。第二步,将求出的多项式特解作为基函数,对三维Burgers方程中剩余的非线性项进行迭代求解。与时空多项式函数作为基函数对三维Burgers方程进行直接求解相比,该算法简单易行,得到的近似解精度非常高,算法极其稳定,对于教学过程中提高学生的编程能力,加深对高维Burgers方程的理解能力以及Burgers方程的实际应用具有重要意义。     

基于简化的应变梯度理论下Kirchhoff板模型边值问题的提法及其应用北大核心CSCD

考虑应变梯度和速度梯度的影响,建立薄板控制微分方程及给出其边值问题的提法,修正了前人给出的薄板角点条件.采用Levy法,给出受分布力作用下简支板的挠度及自由振动频率的解析解.通过与文献中分子动力学数据对比,验证了该文模型的有效性并提出校核材料参数的一种方法.研究结果表明,增大弹性地基和应变梯度参数可以有效提高板的等效刚度,而速度梯度参数则相反.该文提出的板的边值问题为研究薄板在复杂支撑边界及外荷载等条件提供了理论依据.同时,有望为其有限元法、有限差分法和基于能量原理的Galerkin法等数值方法提供理论依据.     

Elastic Analysis of Anisotropic Functionally Graded Rotating Disks With Non-Uniform Thicknesses北大核心CSCD

研究了任意梯度变化的变厚度各向异性转动圆盘的弹性问题.假设圆盘绕刚性轴匀速转动,其材料性能和厚度沿径向任意梯度变化.考虑圆盘在中心转轴处受位移约束,外侧自由,根据各向异性转动圆盘的平衡微分方程,得到关于径向应力的Fredholm积分方程,继而通过对Fredholm积分方程进行数值求解,得到结构的位移场和应力场.对具体梯度变化情况仅需代入相应梯度变化进行求解即可.数值算例部分,通过假设厚度、弹性模量等参数为特殊的幂函数形式,将由Fredholm积分方程求出的数值解与对应的精确解进行对比,以及针对常见的Voigt模型,将由该方法算得的数值解和ANSYS有限元计算结果进行对比,验证了该方法的准确性和精度.其次,针对Voigt模型,重点分析了厚度变化、材料性能梯度参数、各向异性度等对应力场和位移场的影响.提出了针对材料性能和厚度沿径向呈任意梯度变化的圆盘结构弹性分析方法,将为优化功能梯度圆盘的结构和材料参数、有效调整构件应力分布、提高结构安全性,提供强有力的工具;算例分析结果对功能梯度圆盘在复杂条件下的结构安全设计有重要的理论指导意义.     

Impact Responses of Prismatic Lithium-Ion Battery Based on the Membrane Factor Method北大核心CSCD

Aimed at the internal short circuit problem due to large deformation of the prismatic lithium-ion battery cell under impact loadings, a simplified battery model was first established. Then the motion equations of velocity and displacement based on the membrane factor method were proposed. With the effects of the face-sheet thickness and the densification region on the normalized final deflection, impact response characteristics of prismatic battery cells were investigated in detail. The results show that, the improved motion equations involving the membrane factor can reflect the dynamic response mechanisms of the prismatic battery cell under impact loadings, and the large deflection under high-speed impact can be predicted. With the increase of the face-sheet thickness, the deflection of the battery cell’s lower part decreases obviously. However, the densification region expands with the face-sheet thickness. The deflection and the densification region of the cell’s lower part both increase with the inner core density of the battery. This proposed impact model provides a theoretical guidance for the multi-functional integrated dynamic design of prismatic battery cells. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Schrödinger方程的连续时空有限元方法

本文讨论Schr?dinger方程的连续时空有限元方法,通过引入相应的时空投影算子,利用实部虚部分离技巧,得到了变量u在时间节点处的L²范数,以及u和u_t的全局L²(H¹)和L²(L²)范数意义下的最优误差估计结果.该文的结论对进一步探索和设计Schr?dinger方程的数值算法是有益的.